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Abstract 

To  successfully  tune  the  DARHT  II  transport  beamline 
requires  the  close  coupling  of  a  model  of  the  beam 
transport  and  the  measurement  of  the  beam  observables  as 
the  beam  conditions  and  magnet  settings  are  varied.  For 
the  ETA  II  experiment  using  the  DARHT  II  beamline 
components  this  was  achieved  using  the  SUICIDE 
(Simple  User  Interface  Connecting  to  an  Integrated  Data 
Environment)  data  analysis  environment  and  the  FITS 
(Fully  Integrated  Transport  Simulation)  model.  The 
SUICIDE  environment  has  direct  access  to  the 
experimental  beam  transport  data  at  acquisition  and  the 
FITS  predictions  of  the  transport  for  immediate 
comparison.  The  FITS  model  is  coupled  into  the  control 
system  where  it  can  read  magnet  current  settings  for  real 
time  modeling.  We  find  this  integrated  coupling  is 
essential  for  model  verification  and  the  successful 
development  of  a  tuning  aid  for  the  efficient  convergence 
on  a  useable  tune.  We  show  the  real  time  comparisons  of 
simulation  and  experiment  and  explore  the  successes  and 
limitations  of  this  close  coupled  approach. 

I.  INTRODUCTION 

Particle  accelerators  use  magnetic  transport  to 
propagate  the  electron  beam  from  the  accelerator  to  the 
target.  In  the  case  of  the  DARHT  II  accelerator  a  discrete 
set  of  solenoids  and  quadrupoles  performs  this  transport. 
In  the  design  phase  of  the  project  magnet  locations, 
profiles,  and  field  strengths  are  determined  based  on 
various  design  transport  codes  [1].  However,  when  the 
accelerator  becomes  operational  the  adjustment  or  tuning 
of  this  transport  line  will  be  required. 

For  high  rep  rate,  low  energy  accelerators  where  many 
reproducible  pulses  are  available  and  there  are  sufficient 
diagnostics  the  tuning  process  can  be  accomplished  by  the 
operator  using  his/her  understanding  of  transport  and  the 
feedback  received  from  each  shot.  The  Experimental  Test 
Accelerator  II  (ETA  II)  at  the  Lawrence  Livermore 
National  Laboratory  with  a  pulse  rate  of  1  Hz,  energy  of  5 
MeV,  current  of  ~2kA,  and  pulse  duration  of  ~50ns  meets 
the  majority  of  these  criteria  although  diagnostics  are 
relatively  sparse.  The  DARHT  II  accelerator  with  a  pulse 


rate  in  the  neighborhood  of  3  pulses/h,  energy  of  ~18 
MeV,  current  ~2kA,  and  pulse  duration  of  2psec, 
however,  cannot  be  tuned  in  this  manner. 

A  more  practical  approach  for  transport  tuning  in  this 
environment  is  to  employ  a  transport  model  closely 
coupled  to  the  transport  as  a  tool  for  operator  guidance. 
In  the  recent  experiments  with  the  DARHT  II  transport 
beamline  deployed  on  the  ETA  II  accelerator  [2]  this 
approach  was  attempted  using  the  LAMDA  code  [3]  and 
the  FITS  code  [4]  for  the  modeling.  This  paper  reports 
the  successes  and  failures  of  the  FITS  modeling 
experience. 


II.  FITS  Model  of  the  DARHT  beamline  in 
SUICIDE 

The  portion  of  the  DARHT  beamline  to  be  modeled 
starts  at  an  aperture  plate  which  defines  the  beam  radius 
and  limits  the  current  transmitted.  After  this  there  are  two 
solenoids,  five  quadrupoles,  and  then  three  further 
solenoids  leading  to  the  target  as  illustrated  in  Fig.  1. 
Insertible  foils  on  the  sketch  show  where  optical  data  on 
beam  size  is  available. 

The  FITS  (Fully  Integrated  Transport  Simulation) 
model  is  a  beam  transport  model  with  a  graphical  user 
interface  and  copious  diagnostics  which  runs  under  the 
SUICIDE  (Simple  User  Interface  Connecting  to  an 
Integrated  Data  Environment)  system.  The  beam 
dynamics  model  is  selectable  from  several  options 
including  envelope  models  and  particle  models  with 
various  dynamics  effects  included  or  excluded.  For  the 
transport  of  the  beam  through  quadrupoles  the  particle 
models  are  used  exclusively. 


ETA 

Aperture  Foil 

IDIDID 


Quadrupoles  Solenoids 


CQV  CQX 
CQH  CQW 


Figure  1.  Sketch  of  the  DARHT  II  transport  beamline. 


The  FITS  model  provides  for  the  calculation  of  beam 
transport  for  a  single  configuration  and  for  the  calculation 
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of  transport  with  parameter  values  swept  through  a 
prescribed  space.  It  is  capable  of  searching  through  a 
multiple  parameter  space  to  match  observed  beam  radii 
with  FITS  predictions.  This  capability  is  used  to 
determine  beam  injection  conditions. 

The  FITS  model  is  closely  coupled  to  the  ETA  II 
accelerator  in  that  it  reads  the  magnet  current  directly 
from  the  control  system.  The  user  can  grab  the  magnet 
tune  and  then  vary  one  magnet  to  get  a  prediction  of  the 
consequences  of  this  variation. 

FITS  has  been  used  extensively  to  determine  the  beam 
parameters  at  injection.  An  early  FITS  triumph  was  the 
demonstration  that  the  beam  energy  was  ~5.0MeV  and 
not  the  5.8  MeV  previously  believed.  Since  this  was  a 
model  prediction  it  required  a  year  of  digestion  and  two 
other  energy  measurements  to  convince  experimentalists 
that  indeed  FITS  was  correct. 

We  have  found  that  for  transport  through  a  single 
magnetic  element  FITS  agreement  with  experiment  can  be 
made  excellent.  Flowever,  when  multiple  elements  are 
involved  the  accumulation  of  uncertainties  renders  the 
predictions  problematic. 

Ill,  Beam  Size  Data 

The  primary  connection  between  theory  and  experiment 
is  the  beam  size.  There  are  many  code  diagnostics  used  to 
parameterize  beam  radius  and  diameter.  For  envelope  and 
particle  modeling  the  RMS  radius  is  often  preferred  while 
other  transport  codes  utilize  an  edge  radius. 

The  observed  beam  size  is  usually  detemiined  from 
analysis  of  light  generated  by  the  beam  striking  an 
insertible  probe  in  the  beamline.  These  images  are 
characterized  by  several  uncertainties,  the  most  significant 
is  the  optical  background  at  large  radii.  Diagnostics 
which  are  sensitive  to  this  background  level  include  the 
RMS  radius  (the  favorite  of  envelope  modelers)  and  the 
MTF  or  Modulation  Transfer  Function  radius  (favored  by 
radiographers).  Varying  assumed  background  levels  can 
change  the  calculated  radius  by  up  to  ~50%. 

For  these  experiments  the  beam  RMS  radius  is 
calculated  from  the  digitally  captured  beam  image  by 
fitting  a  rotated  two  dimensional  Gaussian  to  the 
smoothed,  cropped  image.  This  gives  a  major  radius, 
minor  radius,  and  angle  of  rotation.  The  fit  automatically 
determines  a  background  or  base  level  for  the  image. 

Although  this  is  a  well  defined  algorithm  which 
produces  reasonable  results  it  must  be  noted  that  other 
reasonable  interpretations  are  possible.  The  beam  radius 
measurement  for  large  beams  is  further  complicated  by 
the  finite  size  of  the  insertible  probe  and  the  presence  of 
return  current  bars  which  shadow  the  image.  For  small 
beams  the  light  generation  mechanism  on  the  probe  is  in 
question  leading  to  further  ambiguity.  In  spite  of  these 
concerns  the  RMS  radius  as  determined  above  will  be 
used  for  comparison  with  the  models. 


IV.  BEAM  ENTRANCE  CONDITIONS 

The  modeling  of  beam  transport  requires  that  the  beam 
entrance  conditions  to  the  transport  section  and  the 
mechanicaFmagnetic  configuration  of  the  transport 
section  be  well  characterized. 

Beam  entrance  conditions  include  the  macroscopic 
parameters  such  as  the  current,  I;  the  radius,  R;  the  rate  of 
change  of  the  radius,  R’;  the  emittance,  s;  and  the  energy, 
E.  The  actual  microscopic  distribution  of  particles  in 
phase  space  may  be  important  but  is  not  measured  - 
simple  Gaussian  or  flat  distributions  are  assumed. 

For  the  DARHT  transport  on  ETA  the  entrance  point  is 
defined  at  an  aperture  where  the  beam  current  is  reduced 
from  ~1800A  to  ~600  amps.  This  aperture  has  an 
opening  of  0.75”  diameter  which  for  a  flat  profile 
translates  to  Rrms  =  0.6735  cm.  The  beam  current  I  is 
determined  by  the  next  beambug  after  the  aperture.  The 
beam  energy  E  is  based  on  energy  analyzer  measurement 
rescaled  to  account  for  the  FITS  observations  of  the  actual 
lower  energy  than  implied  by  the  historical  calibration  of 
the  analyzer. 
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Figure  2.  Determination  of  beam  entrance  conditions  by 
FITS  modeling  of  an  EF3  magnet  scan.  Measured  size 
shown  as  boxes  with  FITS  calculation  shown  as  a  lines. 


The  remaining  beam  entrance  conditions  of  R’  and  e  are 
determined  by  scanning  the  magnetic  field  of  the  next 
solenoid  and  observing  the  beam  size  on  the  next 
available  insertible  foil.  The  FITS  mn  showing  the 
agreement  between  prediction  and  data  is  shown  in  Fig.  2. 
The  large  R’  was  of  concern  because  the  beam  was  tuned 
to  a  minimum  size  at  the  aperture  and  an  R’  of  0  was 
anticipated.  A  FITS  scan  of  entrance  conditions 
addressing  the  aperture  demonstrated  that  the  minimum 
beam  size  on  aperture  was  not  a  waist  but  rather  a 
significantly  diverging  beam,  see  Figure  3. 


minimum  at  apcrtun- 

Figure  3.  Entrance  conditions  showing  how  beam  is  at  a 
minimum  but  not  at  a  waist. 
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V,  Transport  and  Predictions  -  Few  Magnets 

As  the  beam  was  being  brought  online  a  number  of 
prediction  of  magnet  currents  were  made  to  assist  the 
operator  in  setting  individual  magnets.  Transport  through 
the  first  two  solenoids  was  well  predicted.  A  tune  was 
developed  in  FITS  to  transport  the  beam  through  the 
quadrupoles.  Indeed  the  transport  of  a  round  beam  was 
observed  at  about  this  setting.  We  were  able  to  predict 
tunes  reasonably  well  which  produced  a  round  beam  at 
the  end  of  this  transport.  When  solenoid  magnet  currents 
were  varied  we  could  predict  the  rate  of  rotation  induced 
in  an  oblong  beam  from  the  quadrupoles.  The  closely 
coupled  FITS  model  aided  in  the  initial  transport  of  the 
beam  down  the  pipe. 

A  particularly  impressive  success  came  when  we  were 
unable  to  predict  the  beam  size  versus  S4  magnet  current. 
The  focusing  occurred  at  a  current  approximately  50% 
above  that  predicted  by  FITS.  After  considerable 
speculation  we  examined  the  magnet  and  found  that  many 
of  the  multiple  power/water  connections  on  the  magnet 
were  shorted  together  thereby  reducing  the  field  to  2/3  of 
the  predicted  value.  The  magnet  was  fixed  at  that  time 
and  thereafter  agreed  with  FITS  predictions. 

VI,  Transport  and  Predictions  -  Multiple 

Magnets 

Attempts  to  model  in  detail  the  full  transport  through 
the  magnet  system  (including  solenoids  and  quadmpoles 
only;  excluding  steering,  kicking,  and  sextupoles)  met 
with  far  more  limited  success.  Due  to  the  sparseness  of 
diagnostics  often  single  component  effects  could  not  be 
isolated.  With  only  one  magnetic  element  energized  the 
beam  at  the  next  diagnostic  port  would  be  too  large  to 
measure  properly.  With  multiple  elements  energized  the 
models  ability  to  achieve  quantitative  agreement  was 
limited. 

The  problem  is  illustrated  with  one  scan  of  a 
quadrupole  magnet  in  Fig.  4  and  Fig.  5.  When 
quadrupoles  are  varied  the  beam  varies  from  oblong  to 
round  to  oblong  at  90  degrees  rotation  as  seen  in  Fig.  4. 
This  behavior  was  indeed  observed  but  the  currents  at 
which  the  specific  profiles  (e.g.  circular)  were  observed 
often  did  not  agree  with  the  predictions.  Moreover  the 
beam  sizes  at  these  significant  points  were  quite  different 
from  predictions  as  seen  in  Fig.  5. 

This  inadequacy  of  the  model  to  predict 
multicomponent  transport  requires  an  investigation  of  the 
sources  of  error. 


Figure  4.  Scan  of  the  CQX  quadmpole  showing  the 
characteristic  behavior  of  ellipse  ->  circular  ->  ellipse. 
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Figure  5.  Comparison  of  CQX  scan  and  FITS  prediction 
showing  agreement  in  nature  but  markedly  different  major 
and  minor  radii  at  each  current. 

VII,  Limitation  of  Transport  Predictions 

A.  Magnets 

A  primary  source  of  error  is  the  magnetic  field  of  the 
the  transport  magnets.  The  ETA  II  control  system  allows 
the  operator  to  set  a  Request  value  for  the  magnet  current 
and  then  provides  a  Readback  value  of  the  measured 
magnet  current.  Knowing  the  Request/Readback  values 
allows  the  operator  to  return  to  a  particular  recorded  tune. 
Unfortunately;  these  values  do  not  accurately  represent 
the  currents  in  the  magnets.  In  specific  instances  either 
the  Readback  or  the  Request  value  was  close  to  the  actual 
value.  In  other  cases  neither  value  was  correct.  Late  in 
the  experiment  linear  least  square  fits  were  performed  to 
remove  this  ambiguity.  Unfortunately  the  errors  in 
current  measurements  were  found  to  be  non-repeatable. 
The  notion  of  correcting  bad  diagnostics  by  calibrations 
after  the  experiments  is  fundamentally  flawed;  these 
problems  must  be  addressed  prior  to  operations.  This 
problem  potentially  introduces  error  in  the  ~5%  level. 

The  quadrupole  magnets  contained  iron  and  therefore 
had  residual  fields  which  were  not  measured  and  therefore 
not  compensated  for.  Prescriptions  were  put  forth  for 
minimizing  these  fields  but  with  insufficient  current 
available  on  the  power  supplies  this  could  not  be 
implemented.  The  errors  introduced  are  estimated  as 
~5%. 

The  magnet  profiles  were  based  on  model  calculations 
performed  prior  to  construction  of  the  magnets. 
Subsequently  the  model  and  measured  profiles  have  been 
compared.  While  they  are  quite  close  errors  in  the  ~1% 
range  can  be  present  -  particularly  in  the  case  of  ETA  II 
where  the  magnets  were  being  operated  at  lower  than 
anticipated  currents  due  to  the  lower  beam  energy. 

Finally  magnet  installation  was  not  always  in  accord 
with  the  initial  design.  The  S4  was  installed  without  the 
endcap  due  to  constraints  on  the  pipe  size.  This  field 
distortion  is  calculated  to  be  <~1%. 

In  short  there  are  multiple  small  problems  with  the 
magnetic  field  calculation  (even  without  higher  order 
terms)  which  compromise  the  ability  to  predict  magnetic 
transport. 
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B.  Beam 

The  beam  entrance  conditions  are  determined  from 
observations  of  beam  propagation  with  the  first  solenoid 
magnet  scanned.  These  conditions  are  then  used  for  the 
experiment  on  subsequent  days.  The  presumption  of 
constant  entrance  conditions  is  not  justified  for  the  ETA  II 
operating  conditions,  variations  of  upwards  of  ~5%  may 
be  present. 

Only  the  macroscopic  beam  parameters  are  measured 
with  the  FITS  fitting  technique.  However,  microscopic 
details  such  as  the  true  beam  profile,  energy  distribution, 
and  velocity  distribution  may  influence  the  transport. 
Various  distributions  within  FITS  have  been  run  and  for 
most  cases  remarkable  agreement  between  models  is 
observed  but  there  is  always  some  slight  divergence;  say 
~2%.  This  is  especially  true  for  an  expanding  beam  such 
as  is  being  launched  for  these  experiments. 

Oftentimes  one  can  obtain  agreement  between  the 
model  and  observation  by  arbitrarily  adjusting  entrance 
conditions.  While  this  is  at  times  gratifying  it  constitutes  a 
postdiction  and  not  a  prediction  of  beam  behavior. 

C.  Model 

The  beam  dynamics  models  in  the  FITS  code  have 
very  real  limitations.  The  envelope  model  is  only  good 
for  m=0  propagation  and  cannot  transport  a  beam  through 
quadrupoles.  The  particle  models  have  a  primitive 
implementation  of  space  charge.  All  models  are  “single 
disk”  and  ignore  the  head  to  tail  beam  and  field  evolution. 

To  validate  or  vilify  the  models  one  must  pursue 
comparisons  with  other  models  as  well  as  data.  A  detailed 
comparison  with  the  LAMDA  code  is  being  pursued  - 
preliminary  comparisons  showed  agreement  in  principle 
but  detailed  comparisons  are  required. 

D.  Diagnostics 

As  noted  earlier  the  diagnostics  for  determining  beam 
size  have  significant  limitations  in  particular  for  large  or 
small  beams  or  images  with  significant  background  levels 
(that  is  -  all  images).  A  key  parameter  for  comparison  of 
code  and  measurement  for  quadrupole  scans  is  the  current 
at  which  the  beam  is  round  -  the  determination  of 
roundness  depends  on  the  calibration  of  the  orientation  of 
the  insertible  probe  which  varies  with  each  insertion. 
Beam  orientation  errors  range  around  5%. 

VIII,  Conclusions 

Multiple  problems  with  these  experiments  have  been 
identified  and  solutions  can  be  proposed.  However  there 
is  a  fundamental  question  of  the  “predictability”  of  beam 
transport  which  must  be  asked.  A  large  number  of  small 
uncertainties  in  the  dynamics  of  beam  transport  have  been 
identified.  One  can  systematically  try  to  reduce  these 
uncertainties.  However;  there  may  be  fundamental 
limitations  on  the  predictability  of  the  complex  system. 
This  problem  is  reminiscent  of  weather  prediction  where 


the  forecast  for  one  day  (one  magnet)  is  pretty  good  but 
the  seven  day  (seven  magnet)  prediction  is  not. 

One  tack  is  to  try  to  reduce  the  uncertainties  highlighted 
in  section  VII  for  the  four  areas  mentioned.  After 
completion  of  this  series  of  experiments  we  have  begun  to 
explore  this  course  with  an  improved  optical  diagnostic 
probe.  If  we  can  model  the  transport  we  can  tune  with  a 
limited  number  of  diagnostic  ports. 

Another  view  is  that  the  model  will  give  qualitative 
guidance  to  the  tuner  but  not  a  quantitative  predictive 
capability.  This  approach  starts  with  a  “predicted”  tune 
but  then  provides  the  operator  primarily  with  guidance  as 
to  which  knob  to  turn  and  in  what  direction  rather  then 
giving  a  definitive  setpoint.  This  approach  will  require 
more  shots  on  target  than  a  fully  calibrated  model  would 
necessitate. 

Present  efforts  are  focused  on  determining  where  the 
predictive  shortfalls  occur  and  to  what  extent  they  can  be 
remedied.  We  will  continue  to  research  ways  to  improve 
the  predictive  capability.  If  the  opportunity  arises  we  will 
perform  additional  experiments;  particularly  if  better 
magnet  control  /  readout  can  be  achieved.  Simultaneously 
we  are  reviewing  the  heuristic  algorithms  to  be  employed 
for  optimal  convergence  on  appropriate  DARHT  II 
transport  tunes. 
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